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Disordered RL — C networks are known to be an adequate 
model for describing fluctuations of electric fields in a random 
metal-dielectric composite. We show that under appropriate 
conditions the statistical properties of such a system can be 
studied in the framework of the Efetov's non-linear a— model. 
This fact provides a direct link to the theory of Anderson 
localization. 

PACS numbers: 73.20.Fz, 73.20Mf, 05.45.Mt 



Optical properties of random metal-dielectric films 
(also known as cermets or semicontinous metal films) at- 
tracted a lot of research interest recently, both theoreti- 
cally and experimentally, see W-0| and references therein. 
It was discovered that for metal concentrations close to 
the percolation threshold the absorption of microwave ra- 
diation in such materials fluctuates anomalously. In turn, 
these anomalous properties were traced back to high lo- 
cal field fluctuations detected in such compounds. A very 
insightful approach to the problem yj-IJ] is to represent 
the system as a large random network made of capaci- 
tors C and inductances L, the latter being in series with 
a weak resistor R. The network description naturally 
arises when discretizing the equations satisfied by the 
scalar potential of the electric field. The capacitors here 
are to model dielectric bridges, whereas isolated metallic 
granulas are indeed characterized by almost purely in- 
ductive response for frequencies w of radiation such that 



W-, 



< w ^ 



U)r, 



with ujp being the plasma frequency and 



lJt being the plasmon relaxation rate ||l|^ . 

For frequencies close to Uq = 1/(LC)^/^ an electromag- 
netic response of such a network is dominated by reso- 
nance effects as long as losses are small, i.e. the quality 
factor Q = {L/Cy^'^R~^ is large. The resonance fre- 
quencies can be determined as (generalized) eigenvalues 
of some linear lattice operator arising when solving the 
system of Kirchhoff equations J2j '^ij {^i ~ ^i) = for 
on-site potentials Vi = u(ri) [0jH|. Here (Jij is the con- 
ductance between pair of nodes r^ and Vj if two nodes 
are connected by a direct bond and Uij = otherwise. 
In the simplest case, one can think of the network being 
connected to AC voltage by two external leads attached 
to lattice nodes with the coordinates r^ and r^ , the cor- 
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responding potentials being va — e~*"* and vb — 0, re- 
spectively. Omitting the common time-dependent factor 
it is easy to see that the amplitude of the potential v{Yi) 
at an internal lattice node r^ is given by: 



-(ri) = E(^"0./^^ 



(1) 



D, 




(l-%)c 



In a random RL—C network each nonzero conductance 
Uij at frequency / = uj /ii: is equal to either uq — iCoj or 
(Ti = {R + iLuj)^^ , with a specified probabilities (in what 
follows we concentrate on the case of equal probability 
for finding L and C bonds in the network). Then it is 
convenient to introduce "symmetric" variables ft.;,- such 
— 1 if Uij ~ (To and hij = 1 if CTj 



that hi 
that 



'jj 



= eijh ([o-Q + CTi] + [cti - CTo]%), with 



: (Tl, so 

e,y = 1 



for directly connected nodes and e^ = otherwise. In 
terms of these variables we can write D = H—XW, where 



W,, 



{Z + CAt + eBi)Sij - (1 - Sij)ei 



(2) 



Hji — \ h. 



hsi + 



kyii 



Sij - (1 - Sij)h^j (3) 



with Z — "^^ Cij standing for the coordination number 

of the lattice and hij — hijCij. The frequency-dependent 
parameter A is defined as 



A = 



(To - CTi 



S-O-i^-^-i <^' 



where we have made use of a; « wq and Q ^ \. 

We see, that statistics of the scalar potential v{t) (and 
hence of the electric field £ij proportional to the voltage 
difference v{ri) — v{vj) on the bond ij) is determined ba- 
sically by properties of the operator H . The operators 
of such a type acting on a lattice were suggested to be 
called Kirchhof Hamiltonians (KH) in [Q . Off-diagonal 
entries of such a Hamiltonian assume random values ±1 
for directly connected nodes. This property makes KH 
to be , in a sense, similar to a tight-binding Hamiltonian 
describing the motion of a quantum particle on a disor- 
dered lattice with an off-diagonal disorder. The latter 
model is a paradigmatic one in the theory of Anderson 
localization. That kind of analogy first discussed in R] 
led the authors to relating the anomalous fluctuations of 



electric fields to localized properties of the correspond- 
ing eigenfunctions. Further numerical and experimental 
work confirmed the qualitative validity of the suggested 
picture. 

At the same time, the question to which extent one can 
push forward the analogy between the Anderson model 
and the Kirchhof Hamiltonian is far from being trivial. 
Indeed, the KH has a specific feature: the diagonal en- 
tries Hii are strongly correlated with the off-diagonal 
ones Hi^j . It is known that correlations of various kinds 
can substantialUy modify the localisation behaviour, see 
e.g. |6|. Therefore, it is highly desirable to find an ade- 
quate approach allowing to shed more light on the ques- 
tion of equivalence between the models. 

The main goal of the paper is to show that the equiv- 
alence indeed exists and the unifying concept is pro- 
vided by the so-called Efetov's supermatrix non-linear 
cr- model (ENSM) 0. The latter model is known to be 
the most powerful tool in understanding the fluctuation 
phenomena in disordered conductors last decade, see e.g. 
g. As a matter of fact, we derive ENSM from a version 
of the Kirchhof Hamiltonian and thus provide a regular 
analytical background for the quantitative description of 
statistical properties of the semicontinous films. 

To derive ENSM from a microscopic random Hamil- 
tonian one has to exploit some large parameter which 
physically controls the strength of the disorder. The ex- 
perience of dealing with the usual tight binding models 
suggests that a role of such a parameter can be played, 
e.g. by a large radius of connectivity b (i.e. the large co- 
ordination number Z '^ b'^) [plfO|. Formally, we consider 
a d— dimensional lattice of a linear size L with unit lattice 
spacing and a connectivity radius b. To facilitate book- 
keeping of terms of the different order it is convenient to 
redefine e^ -^ e^Z'^l"^ where e^j = 1 for jr^ — rj| <b 
and eij — otherwise. The radius of connectivity b is 
chosen to satisfy 1 <^ b <^ L. Both inequalities are im- 
portant: 5^1 allows one to map the problem to ENSM 
, whereas 6 <C i is necessary to ensure the adequate de- 
scription of effects of the Anderson localization. Indeed, 
as is shown recently g a full-connectivity LC— network 
with b = L can be mapped on the zero-dimensional ver- 
sion of the ENSM, which precludes the localization effects 
to be taken into account. 

To demonstrate the mapping it is instructive to ad- 
dress the simplest nontrivial correlation function of the 
potentials: Ci(J7,r) = (w*j(ri)ui^2(ri)) where we intro- 
duce the frequency difference fl oc (wi — ll)2)/'2^o "^ 1 
and the brackets stand for the disorder averaging |ll[ . 

Our starting expression is: 



Ci(ri,r) = ( 2J ^AkiSAkz [JiAki - K) {hAk2 
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(5) 



To perform the disorder average we follow the stan- 
dard procedure and represent the matrix element of the 
resolvent {H — XW)^"^ in terms of the Gaussian integral: 



1 



H - (ReA ± iT/2)W 



±i 



ik 



N 



\{d^i{±) 



.1=1 



s:i±>k{±) 



N 



exp{±- J2 *™(±) [WmniReX ± lT/2) - Hmn] *„(±)} (6) 



over 4-component supervectors ^;(±), 






,Si{±) 



Si{±) 
Xi{±)\ ^^ dsidsl 






(7) 



dx*idxi 



with components s;(-f), Si(— ); I = 1,2,...,A^ being 
complex commuting variables and Xi(+)jXi(~) forming 
the corresponding Grassmannian parts of the supervec- 
tors *;(±). 

To facilitate the presentation, it is appropriate to an- 
ticipate few facts whose validity can be verified by the 
same method as presented below. First of all, after av- 
eraging the double sum in the expression Eq. (m is dom- 
inated in the limit 1 <C & "C i by the diagonal terms 
with indices fci — fc2- Another fact which is useful to 
exploit from the very beginning is that all the resonance 
frequencies are concentrated in an interval of the order of 
5lo /uoq ^ Z^^'"^ around uj — wq, so that the typical spac- 
ing A between the neighbouring resonances is of the order 
of A -- u}q/{NZ'^/'^), with N ^ L"^ being the total num- 
ber of resonance frequencies. We anticipate nontrivial 
correlations occuring on the frequency scale wi — cj2 ^ A 
11 . For this reason we scale ReAi,2 = (r ± n/2N) /Z^/^ 
considering both r and f2 to be of the order of unity. By 
the same reasoning we consider the losses to be small 
enough to ensure that 7 = T{NZ'^/'^) is of the order of 
unity. Physically this requirement means that a typical 
resonance width ujqT is considered to be comparable with 
the typical resonance spacing A. 

With these facts in mind, we can easily average the 
products of the resolvents over hij — ±1/Z^' ■^e^ in the 
limit A^ 3> ^ 3> 1. All further steps follow the method 
used in g (cf. ||l^,|l2[ ) adopted to the present model. We 
have: 






Yld^i 



^K-M+Ki+>kH (8) 



^""p 1 1^(" + '^^ X! ^™" (*™ ~ *») (*'» " *« 



X exp ■ 



ik2 J 



rii.n ) 



Here the integration goes over the 8— component super- 
vectors $ 



I -(*!(+), vl/t(_)) and 



/C($a,$6) = 



(9) 



-Air I 



When deriving Eq.(|l5|) we have used Eq.([lO|). 

Given the form of the kernel Eq.(H) and exploiting 



[K - H) A i^a - $h) + [($1 - H) A ($, - $b)J ZE„[e-n™^ = 1 one can find a space-independent so- 



with A = diag{l,f,l, 1,-1, -1,-1,-1). Only those 
terms are left in the exponent which will later on con- 
tribute to the final expressions in the discussed limit. 

Next step is to use the following functional Hubbard- 
Stratonovich transformation (cf. [|lO|Jl2| ,|5|] ) : 

exp|-^^e„„/C(<I>,„,$„)| = /'p(g)e*S".iS'"(*™) 

l^ m.n ) 

^^"""^ 1 16^ [^"^]™" / «'*«^*65m($a)C($a,$b)g«(*b) 
l^ nin 

where [e^^l is the matrix inverse to the matrix e — 
[cij] |ij=i,....,Ar and the kernel C(<I>a,$b) is , in a sense, 
the inverse of a (symmetric) kernel IC{^a, ^b)- 



I 



(i$/C($a,$)C($,$6)=(5($a,$b), 



(10) 



with 6{^a, ^b) playing the role of a S— functional kernel 
in a space spanned by the functions g(<&). 

With the help of these relations one easily brings each 
term of the sum in Eq.(||) to the form: 

J V{g)T+,_ [9^]:F.,+ [gk] cxp [C{g} + dCi{g}] (11) 



where 



^±Mfj] 



/d$.s*(±)s(T)et9(*) 



■Y^[e^^]mn / d$ad*b5m(*a)C($a, $6)5n($6) , 

.{n + i-f) 



(12) 

(13) 



Z 
16 



SCi{g} 



ANZ 



E 



(14) 



tnKn 



{<^l - $t^) (<;>,, - $, 



/ d$„d$„ exp I [g„^{<^rn) + gni'^n)] 

and we restricted ourselves by the leading order term in 
5Ci{g} which is only a small correction to C{g}. 

Next step is to evaluate the functional integral over 
g{^) by the saddle-point method, justified by two large 
parameters: Z and N . The saddle-point configuration 
gm (^) can be found by requiring the vanishing variation 
of the "action" C{g} and satisfies the following system 
of equations: 



lution gn\'^) = 3^''^$) to equation Eq.(|l|) : 

g^'H'^a) = 4(r - Gi)($lA$a) + 4zG2($l$a) + «($IA$J2 

(16) 

provided the real coefficients Gi, G2 are solutions of the 
system of two conjugate equations: 

G2±iGi^ I dMexpj±^u(r-Gi±iG2)-y| (17) 

jHints to verifying such a solution can be found in |q|. 

For further analysis it is very important that a solution 
to equations Eq. ([l7|) exists for arbitrary — 00 < r < cxo 
such that G2(r) > 0. Actually, the mean density of res- 
onances is merely given by p(r) = -G2(r) [g. 

The most important consequence of the existence of 
the solution ^('^^(^a) in the form Eq.(|l6|) with G2 7^ is 
actually the simultaneous existence of a whole continous 
manifold of the saddle points parametrized as: 



gT{^a)=9^'\f<^a), 



with f'Af = A 



(18) 



If it hadn't been for the condition G2 7^ all these 
solutions would trivially coincide: gri^) = g^'^H^) ^^^ 
any T defined as above. In the actual case presence of 
the combination $a$a which is not invariant with re- 
spect to a transformation <I>q -^ T^a ensures the ex- 
istence of the mentioned manifold. This fact is just 
a particular manifestation of the phenomenon of spon- 
taneous breakdown of symmetry. Different nontrivial 
solutions are actually parametrized by the superma- 
trices T which are elements of a graded coset space 
UOSP{2, 2/A)/UOSP{2/2) ^ UOSP{2/2). 

As a result, the functional integral over gmi^) is dom- 
inated by "Goldstone modes" slowly changing in space 
and parametrized as: gm {^) = g^^H'^m^) , with the 
matrices T,„ which depend on the lattice site index 
m — 1,...,N. Our next step is to determine the effec- 
tive action for these modes that is jC{gm } + SjCi{gm }■ 

First of all we notice that: 

which allows one to perform the following manipulations: 
d<i>„d$,g(G)($JC(<i>„, $fc).gf )($fc) (19) 



i /"d$giG)($)etff"'(*) =i f d<i>g(''\f„f^^<P)ei 



gt^'C*) 



-AGlStrifnf-^Aifnf-^] A 



C, = {npf / VQ{r)Q^^r,)Q^\vA)e 



CiQ) 



(24) 

In conclusion, we managed to express the correlation 
function of the scalar potentials in terms of ENSM. Sim- 
ilar reduction is possible also for other quantities of in- 
terest. As is well-known pfl, explicit evaluation of the 
integrals of the type Eq.(24) crucially depends on the pa- 
rameter g = 2Tr pDL'^^^. For g ^ oo the integral is dom- 
inated by the constant configuration: Q{r) = Qq and the 
result for Eq.(p^) is very simple: Ci = 2i'Kp{r)/{Vl + 17). 
This is the so-called " zero-dimensional" limit correspond- 
ing to the infinite-range connectivity model H. One 
can take into account weak localization effects finding 
1/g corrections in any dimension d, see e.g. [H|. For a 
quasi one-dimensional lattice one can calculate integrals 
Str [f^^kf^k + f^^kfnA) exactly in the limit f7, ?? ^ |||l. One should be able also 

to study singular parts of higher correlation functions. 
, . , 1 v^q; /'t'-Iat' \\ , -These questions are left for further investigations. 

7rp(S2 + ^7)J^ 2^ Str 1^1^ AlmAj (20) ^^^^ author is grateful to J.-M.Luck and V.M.Shalaev 

for their encouraging interest in the work. The present 



= -4(7rp)2Str (f„iAf„f -iAf„ 

Here we first exploited the saddle-point equation Eq. ( p^ ) 
together with Eq.(p^ and then performed a change of 
variables: T„i^ -^ $ which does not effect the measure 
(i$ because of the (pseudo) unitarity of the matrices T. 
Then the integral can be readily evaluated with help of 
the explicit form Eq.([l6|) (see e.g. examples of similar 
calculations in ||5|,[l2|) and brought to the final form by 
employing the cyclic permutation and the mentioned re- 
lation between G2 and the density of resonances p{r). In 
the very same way we also find: 



SCi{g} 



{fl 



rn<n 



The pre-exponential factors J-'-^-_-[gi] are calculated anal- 
ogously and are given by 



.F+ 



[9i\ = T^P 



^AfJ ' ; .F_,+ [5fe]=7rp[f^-iAfJ ' 
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where the indices of supermatrix elements are inherited 
from the structure of the diadic product $ ® <I>'f . 

We see, that calculating the correlation function Eq.(P) 
in the limit L ::^ 6 ^ 1 amounts to evaluating the follow- 
ing integral over the set of supermatrices Qm = T^^^kTm'- 



Cr 



{npf 



^CAk 



N 



WdQi 



.1=1 



(0) 



Z{^P? 



N 

E 



Q^^Q^^^e^(« 



oL^^/7J2^;r: 



(21) 
(22) 



N 



m— 1 



47V 



The action Eq.(p2) is actually equivalent to the dis- 
cretized version of the supermatrix nonlinear ct— model 
well-studied in the context of the Anderson localization 
[0,H. Indeed, following |^ notice that the condition 
6 S> 1 ensures slow variation of the matrices Qm with in- 
dex m, so that it is legitimate to pass from the lattice to 
continuum, when the action assumes the standard form: 



C 



(Q) 



np 



drStr 



D 



M' 



2i 



(» + tj) 

V 



Q(r)A 



[10] 
(23)[11] 



where D — irpZ ^ X]|r|<6 ^^ pl^Y^ the role of the effective 
diffusion constant and V = J dr. The large value of 
D (X b^ ensures that a typical spatial scale S, of variation 
of Q(r) is large: ^ ex b^. Thus, for distances of the order 
of b the matrices Q do not change and therefore: 
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